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The short-range defect with reduced symmetry is studied in the framework of kp-approach taking 
into account a matrix structure of potential energy in the equations for envelope functions. The case 
of the narrow-gap semiconductor, with defects which are non-symmetric along the [001], [110], or 
[111] directions, is considered. Resonant state at a single defect is analyzed within the Koster-Slater 
approximation. The bandstructure modification of the alloy, formed by non-symmetric impurities, 
is discussed and a generalized virtual crystal approximation is introduced. 

PACS numbers: 71.20.Nr; 71.55.Eq 

I. INTRODUCTION 

Examinations of electron states, localized at short-range defects, began with the development of the band theory (see 
Refs. in (lj). The approach for the calculation of the photoionization cross-section was formulated in the framework of 
the effective mass approximation with short-range impurity potential (Lukovsky model Q). The short-range acceptor 
states were described with the use of the multi-band approach based on the matrix kp-Hamiltonian and on the scalar 
potential energy U6(r) for the impurity placed at r = ||. The same approach was used Q for the description 
of both the resonant and localized electronic states in the gapless and p-type materials under uniaxial stress. The 
multi-band considerations in jj| and [Q have demonstrated a good agreement with the experimental data for the 
photoionization cross-section of acceptor state and the transformation of impurity states under uniaxial stress (see 
discussions in Q and || respectively). The considerations ^ |[ ^ were restricted by the presumption that defect 
does not mix different bands by itself due to the scalar character of the impurity potential used. To the best of our 
knowledge, a complete description of electronic states at short-range defect, taking into account a matrix structure of 
potential energy, has not been performed yet. This issue is also of interest in connection with the last investigations 
of A3B5 alloys formed by materials with distinctly different parameters (see reviews in Q). 

In the present paper, we evaluate the generalized kp-approach with the matrix potential energy described a short- 
range defect with reduced symmetry. We demonstrate both the appearance of resonant electron state at single defect 
and the significant modification of alloy bandstructure due to the non-diagonal contributions to the matrix potential 
of defect. The consideration is performed in the framework of 6 x 6 Kane model [p|, corresponding to a narrow-gap 
semiconductor, for the cases of defect which is non-symmetric along the [001], [110], or [111] directions. The matrix 
Green's function of the resonant states is obtained in the Koster-Slater approximation and the conditions for a narrow 
resonant peak over the conduction or valence band are discussed. The low-energy electronic states in semiconductor 
alloy are described by the Dyson equation with the self-energy function, written in the Born approximation. Due to 
the weak damping, we arrive to the effective Hamiltonian with the non-diagonal matrix of the extreme energies, i.e. 
we introduce the generalized virtual crystal approximation. As a result, the energy spectrum appears to be anisotropic 
and splitting of the degenerate hole states takes place. In addition, a zero-gap case can be realized with an increase 
of alloy composition. 

The consideration below is organized as follows. In Sec. II we present the theoretical background including the 
zero-radius approach for description of the short-range defect, the Koster-Slater solution for the resonant state at a 
single defect, and the generalized virtual crystal approximation for the low-energy electron states. The results for the 
resonant states and for the alloy bandstructure are discussed in Sec. Ill with the use of the 6x6 Kane model. In 
Sec. IV we make concluding remarks and discuss the assumptions made. 

II. BASIC EQUATIONS 

We consider here the equations for the envelope function and for the single-particle Green's function in the crystal 
with short-range defects. Further, we evaluate the contribution of the resonant states into the density of states and 
describe the low-energy electronic states using the Dyson equation with the second-order self-energy function. 



2 



A. Zero-radius approach 



The general equation for the envelope wave function, ^/ r , which is introduced through the expansion of an exact 
wave function, ^ ^/rUjr, over the Bloch amplitudes of Zth state in the center of Brillouin zone, u lri is written as: 



E 



H u ,H! Vt + / dv'U w {v,v')^i 



(1) 



Here I includes band and spin indexes, Hw is the single-electron kp-Hamiltonian described a motion in the ideal 
crystal. The microscopic potential of defect, 5U r , appears in the kernel: 



U w (r,r') = J drm^Air-r^SUr.Ain-r')^^, 



(2) 



where A(r) = (1/V) ^ k exp(ikr) is the (5-like function; V is the normalization volume and ^ k ... is taken over the 
first Brillouin zone. Here we consider the short-range (localized in the elementary cell) defect placed at r = 0. For 
the case of the low-energy states, with the wavelength higher than the lattice constant, a, we substitute *;' r /=o into 
the integral term of Eq.(l) and this equation takes form 



(3) 



Thus, we have obtained the equations for the envelope wave function with the matrix potential energy Uw A(r) 
determined through the components 



Uw = j dru* lr A(r)5U T u v 



(4) 



moreover Uw = U* n . It is necessary to stress that there are no reasons in order to neglect the non-diagonal components 
of this matrix if the local symmetry of SU r is reduced in comparison to the crystal symmetry. This point is a main 
difference from the standard kp-description of the large-scale potential case and from the previous consideration of 
the deep impurities cited above. 

Using the potential energy J2j A(r — Rj)U f° r the crystal with N im short-range defects in the volume V placed at 
Rj, j = 1 . . . N im , we write the equation for the retarded Green's function in the matrix form 



E + iX - H - A(r - Rj)U 



G E (r,r') = 5(r-T% 



where A — > +0 and 1 is the identity matrix. Next, we transform Eq.(5) into the equation: 

G E (r, r') ~ g E (r, r') + ^ g B (r, R 3 )UG E (R 3 ,r') 



(5) 



(6) 



where the free Green's function, g E (r, r'), determined from: (E + iX — H)g E (r, r') = <5(r — r')l. Thus, the zero-radius 
approach gives us an infinite (if N im — ► oo) set of equations for G E (Rj,r'). 



B. Koster-Slater approximation 



In order to describe a localized (or resonant) state at jth defect, we rewrite Eq.(6) for G E (Rj,r') neglecting the 
neighbor defect contributions: 



l-g E {Rj,Rj)U G E (Rj,^') — 5 E (Rj,r'). 



(7) 



The function g E (Rj,Rj) means the free Green's function with the cut-off singularity at coinciding arguments (i.e. 

— 1/3 

at distances ~ a). This approximation is valid if the radius of states is shorter than n im ; n im — N im /V is the 
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concentration of defects. Substituting G E (Rj,r') determined by Eq.(7) into Eq.(6) we obtain the Green's function 
written in the framework of the Koster-Slater approximation: 

G E (r,r') ~ g E (r, r') + £ g B (r, Kj)U (l - A E tf) _1 g E (Rj , r') (8) 

j 

with A E = g E (r,r). 

The density of states is expressed through G E (r, r') according to the standard formula: p E = 
— Tmtr J drG E (r,r)/(TrV) where tr . . . means sum over diagonal matrix elements. The impurity contribution to the 
density of states, Ap E , is determined by the second term of (8). After the permutation of g E under tr . . ., we use the 
equality J drg E (Rj,r)g E (r,ILj) = — [dg E (Rj,'R,j)/dE} and Ap E is transformed into (see similar formulas in (J, D) : 

Ap E = ^I m tr^u(l-A E uy 1 , A E = LJ2'9 E (P), (9) 

p 

where ^2' p . . . means the summation over the region |p| < fi/a. Here we have also used the free Green's function in 
the momentum representation: g E (p,p') — S p . P 'g E (p) withg E (p) = (E + iX — e — v • p) _1 , where the kp-Hamiltonian 
in p-representation, e + v ■ p, is written through the extreme energy and interband velocity matrices, e and v. It is 
convenient to write g E (p) through the dispersion laws of fcth band, e^p, according to: 

kp ^ 

where the projection operators onto the fcth band, Pkp, are given by the matrix elements: (Pkp)w = ^ip^* ■ 
Here a is the spin index and the column ipp ka ^ is determined through the eigenstate problem: (£+v-p)-0 p fcc ^ = ekpipp?^ ■ 



C. Generalized virtual crystal approximation 

— 1/3 

Consideration of the low-energy electronic states, with the wavelength exceeds n im , is based on the matrix Dyson 
equation for the averaged Green's function. Such equation is obtained from Eq.(6) in the form: 

G E (p)=g E (p)+g E (p)£ E G E (p), (11) 
p 

where the self-energy function, S E , is written in the self-consistent approximation. Using <? E (p) in E E (the Born 
approximation) we write the self-energy function through A E given by Eq. (10) according to E E ~ rii m tj L E \] . 

Since Eq.(ll) is transformed into G E (p) = [gsip)^ 1 — S E ] _1 , we obtain the averaged density of states in the form: 

^=-^5>(^-"- p -^r (i2) 
p 

and matrix £ E reduces the symmetry at p = 0. According to the consideration below (see Eqs.(16)-(20)), the self- 
energy matrix appears to be weakly dependent on E and the damping contributions to £ are also small. Thus, we 
arrive to the effective Hamiltonian: 

=i+t + *-p, £~ n im UAU (13) 

with the non-digonal matrix (e+S) which determines the energies of the band extremes. This Hamiltonian corresponds 
to the generalized virtual crystal approximation for the alloy formed by low-symmetric defects: not only the energies 
of the band extremes are shifted with the alloy composition (see |fl0| ; note, that E oc ni m in the Born approximation), 
but also the symmetry of H p eff) is reduced due to the non-diagonal contributions from S. The characteristic equation 
correspondent to the effective Hamiltonian (13) is written in the form: 

det|e + E + v-p-£| =0 (14) 

and the dispersion laws appear to be anisotropic, with splitted valence bands, like for the case of stressed semicon- 
ductors. 
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III. RESULTS 

Let us turn to Eqs. (9) and (14) for the narrow-gap semiconductor described in the framework of the isotropic 
6x6 Kane model. For the sake of simplicity, we restrict our consideration to the model of defect with the microscopic 
potential, 8U T , which is only non-symmetric along the [001] (case A), [110] (case B), or [111] (case C) directions. 



6x6 Kane model 



The matrices U and A E are determined below using the set of Bloch functions which are expressed according to 
JTl| in terms of the periodic basis \S), \X), \Y), and \Z). Within the above approximations, the potential matrix (4) 
is written in the form: 



U = 



u 



-i 



1 o o 



Ud 



7§ e + 

7§ e - 





ie- 
— ie + 



(15) 



7T 



'\/| ie-_ 
7§ e - 



Here ii c and u v are proportional to the 2 x 2 and 4 x 4 identity matrices: u c = 1(S\8U\S) = lu c and u v = l(X\5U\X) = 
lu v ; the non-diagonal part of U is expressed through the 4x2 matrix, u, where e± = (l±i)/V2. The matricis u A B C 
are proportional to u z = (S\SU\Z) + with (S\SU\X) = (S\SU\Y) = (case A), or u d = (S\SU\X) = (S\SU\Y) ^'o 
with (S\SU\Z) = (case B), or u a = (S\SU\X) = (S\SU\Y) = (S\SU\Z) ^ (case C). 

Since the dispersion laws Sk P are isotropic, we use in Eq.(10) the projection operators which are averaged over the 
angle, so that Pk p are proportional to Sw . Thus, the diagonal matrix A E is determined by the matrix elements 



(A E )u 



1,2 
3-6 



(16) 



where Xk E are weakly dependent on E due to the dominant contribution from |p| ~ h/a. Because of this, we use the 
expansion \k E — — lk(E + e g /2) where the energy —e g /2 corresponds the middle of gap and the coefficients 
and h are given by 



Xr 



1 E'- 



\C a 



V ^ £ a/ 2 + £ iv 



lr 



\C a 



1 E' 



|Cgpl a qp 



1 E 



\c, 



qp I "'qp 



V^(e g /2 + s qp y 



p 

l -Y^(e g /2 + e hp )- 2 . 



V 



(17) 



V 



Here ^2 q . . . means summation over the electron (o) and light-hole (lh-) state while the heavy hole (hh-) state, with 
the dispersion law Eh P , appears in A„ and l v . Note also, that the coefficients Ik have no divergence at |p| — > oo. Eqs. 
(17) are written through the coefficient, a qp , and the normalization factor, C qp , which are determined as follows: 



[Vvf 



3(e qp + s g ) 2 



\C a 



2(V P ) 2 /3 



(18) 



where V is the Kane velocity. The broadening factors in Eq.(16), r]k E , are expressed through the fcth band density of 
states, Pk{E) = (2/V) J2 P S(E - e kp ), according to 



= ? E \ C ip\l qp = E Pi{ E )i Vv E = | y,{\C qp \ 2 a qp ) Eqp=E p q {E) + j Phh (E). 



(19) 
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Thus, the model under consideration involves five phenomenological parameters: the potentials u c . v and u z (case A), 
or Ud (case B), or u Q (case C) in Eq.(15) and the factors A Ci „ determined by Eq.(17). Another values (l c<v in Eq.(17) 
and the functions -q c , VE given by Eq. (19)) are expressed through the parameters of the ideal crystal determined in 
kp-approach. 



B. Resonant states 



We consider here Ap E determined by Eq.(9) taking into account both diagonal and non-diagonal contributions to 
Eq.(15). Below we use 6x6 matrix A E U and we introduce (1 — A B J7) _1 in the matrix form: 



X 



CE^C ^CE^ 



{i-A E uy 



2-C •£ c v 



(20) 



where the 4x2 matrices 



L A,B,C 



are given by the Eq.(15). The components of the reciprocal matrix in (20) are 



determined from the 6x6 matrix equation: 



\ CE u c 



\ CE U 

All R till 



) 


x c 











Using these notations, we write Eq.(9) as follows: 

Ap E = — ^Im [\' CB tr (u c x c + ux vc ) + X' VE tr (u v x v + u + x cv )] . 



(21) 



(22) 



where X' qE = dXn E /dE. Further calculations reduce Eq.(21) to the 2x2 linear equations and result for Ap E takes 
standard form m, \ 



a ^ ^ 2 171 -r 

Ap E = Ira 



(d£ E /dE) K' VE u v 



1 A. VE u v 

where the last term is due to the heavy hole contribution. The function C E is written as 

2^/3, A 

C E — (1 - X CE u c )(l - X VE u v ) - wX CE X VE , w = { 4m^/3, B 

2u% C 



(23) 



(24) 



with the different factors w for the cases A-C under consideration. 

According to Eqs.(16)-(19), C E depends on E weakly and has the small imaginary part. It is convenient to introduce 
the level energy E Q which is given by a root of the linear equation: ReC Eo = 0. Thus, the function C E is transformed 
into C' a {E — E — iT ) where the coefficient, C , and the broadening energy r o , are introduced as follows: 



C = l c u c (l - X v u v ) + l v u v (l - X c u c ) + w(l c X v + l v X c ), 
r o = [r]cE u c + r] VBo u v + w(rj CEo X v + rj VEo X c )]/ C . 



(25) 



Note that rjkm and T a are proportional to the density of states in c- or w-band, if E a > or E Q < —e g . It follows that 
the broadening energy for the w-band resonant state, with E a < —£ g , exceeds r o for the c-band resonant state, with 
E > 0. The same transformation can be performed for the last term of Eq.(23), so that the level, E v , is determined 
from the equation 1 — X VE u v = and the broadening energy is equal to rj VEv /l v . 

If a deep local level appears in the gap, > E QyV > —e g , then the broadening energy is replaced by +0 and the 
contribution to Ap E takes form: 2rii rn 5(E — E D _ V ). The cases E 0>v > or E 0tV < —e g are corresponded to a resonant 
state over c- or w-band and Ap E has the Lorentzian shape with the broadening energy T : 



Ap E 



(E-E )2 + T* 



(26) 



The same result with the level energy E v and the broadening energy rj VEv /l v is valid for the /i/i-contribution. Here and 
in Eq.(24) we have supposed that E D:V or \E Q . V — e g \ exceed the broadening energies for the resonant states over c- or 
w-bands respectively. Both E a and r o in Eq.(26) are determined by the short-range contributions to the factors A c ,„, 
so that neither an absolute value of broadening energy nor a sign of r o , which determines peak- or dip-modification 
in the density of states, are not fixed in the consideration performed. The only dependence on e g (i.e. on hydrostatic 
pressure) can be evaluated from these results within a few fitting parameters. 
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El,2 




£5,6 


A 


n im {u 2 c \ c + u 2 z \ v ) 


—E a + n im (u 2 \ v + UgXc) 


. 2 \ 
Eg T TlimUv 


B 


riim(ui\ c + 4UrfA«/3) 


— Eg + nim(ul\v + Uj,\c/3) 


— Eg + nim{u 2 \ v + U 2 d \ c ) 


C 


n lm (u 2 c \ c + 2u 2 ,\ v ) 


— Eg + riim(ulX v + U 2 \ c ) 


—Eg + n im (ul\ v + u 2 \ c ) 



TABLE I: Components of diagonal matrices -61,3,5 in Eq.(27). 




C. Alloy bandstructure 



Next, we examine modifications of the alloy bandstructure due to contributions into Eq.(14) from the matrix 
self-energy function determined by the Eqs.(15), (20). Using 2x2 matrix notations, we write £ + £ in the form: 



E 1 
i? 3 A 
A+ E 5 



(27) 



where the non-diagonal component, A, is determined for the cases A-C through the characteristic energies S 

A c /\/3 according to: 



n lm u 2 d \ c / V% and <5 



A A = 0, A f 



iS, 



% -(1 + 



(28) 



The diagonal matrices -£1,3, 5 are determined by the components presented in the Table. 

The spin-degenerate dispersion equations are obtained from the Eqs.(14), (27) and (28) as follows: 



(7>Px) 



-(E-E 3 



(E-E 5 ) = (E-Ex) x 



(E - E 3 )(E - E 5 ), (A) 
(E-E 3 )(E-E 5 )-S 2 , (B) 



^(T P ) 2 (E - E 3 ) + ^6T 2 ( PxPy + PxPz + PyPz ) = {E- E x ) [(E - E 3 ) 2 



(C), 



(29) 



where P \ = P 2 + P 2 . According to Eq.(29) the extreme energies, S3, 5 in the case A and E± — E 3 ± \8\ in the case B, 
appear to be linear dependent on the alloy composition [ p~2| . A bowing of u-band extreme energies with increase in 
alloy composition takes place for the case B 



E 3 +E 5 
E± = Z ± 



S3 + S 5 



+ 6 2 



(30) 
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VPJS 7>pje t 



FIG. 2: The same as in Fig.l for the case B. 




'PP. /£ 



FIG. 3: The dispersion laws along the cubic diagonal p||[lll] (pd, solid curves) and p\\OZ (p z , dotted curves) for the cases C 
if £ 3 ,5 = -0.2 eV (a) and £3,5 = -0.3 eV (6). 



due to the non-diagonal contributions in Eq.(27). The anisotropy of energy spectrum takes place if £3 7^ E 5 in the 
cases A and B, or if <5 7^ in the case C. These dispersion laws are shown in Figs. 1-3 for the narrow-gap alloy with 
2m V 2 —25 eV, m is the free electron mass, and with e g =0.25 eV. Here we suppose the negative shift of c-extremum, 
Ei = —0.1 eV, for the all cases A-C and the energies S , 5 are chosen as 0.05 eV. We consider a different order of 
w-bands: E 3 = -0.2 eV, E 5 = -0.3 eV (a) and E 3 = -0.3 eV, E 5 = -0.2 eV (b) for the cases A, B in Figs. 1,2, and 
we suppose £3-6 = —0.3 eV and -E3-6 = —0.2 eV in Figs. 3a and 3b respectively. One can see that the dispersion 
laws appear to be weakly anisotropic in the cases A and B while a visible anysotropy takes place for the case C. We 
note also a substantial modification of the w-band extrema if the splitting energy, \E 3 — E§\, is comparable with e g . 

The zero-gap bandstructure can be realized if E\ = E% or if E\ = E$ for the case A and if E\ = E + for the case B 
and if E\ = E 3 + 5 for the case C. In the vicinity of the cross-point energy, if \E — £a|,'P|p| are less than the energy 
spacing to the lower w-band, the dispersion laws are given by: 



E „ E , + V / (^) 2 /6 + (2/3)(^) 2 , E 1 = E 3 

I ±V P JV2, E x =E 5 [ ' 



F ~. F +T [Pi . 2pj\ E.-Es pl 

E ±P *E 1 ±PJ— —\- + — +— —- (B) (31) 
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E ±p ~ Ei ± —^yP 2 + sign(5){p x p y + p x p z +p y p z ), (C). 

Thus, we have obtained a linear dispersion laws which demonstrate an essential anisotropic behavior in the vicinity 
of the cross-point energy. 

IV. CONCLUDING REMARKS 

In this paper, we have developed the generalized kp-approach for the description of the short-range defects with 
reduced symmetry. We have considered the peculiarities of the density of states due to the resonant state contributions 
as well as those of the alloy bandstructure. In contrast to the previous considerations, here we have used the matrix 
form of the potential energy in kp-equations which leads to the additional interband mixing. In spite of the simplified 
model including a few phenomenological parameters, we have demonstrated an essential modification of the results 
due to the non-diagonal part of matrix potential. 

Next, we discuss the main assumptions used. The above approach is based on the single-particle description of 
the electron states which is valid for many kinds of defects in semiconductors. Although the general equations 

(3)- (12) are based on the presumption of low energy of electrons only, the concrete results are restricted due to the 
specific microscopic potentials used and due to the narrow-gap approximation based on the isotropic 6x6 Kane model. 
The last approximation supposes that far-band contributions do not change the results qualitatively. We have also 
used the low-concentration approximation, supposing that the electronic states at different impurities do not overlap 
in Sec. IIB, and we have used the Born approximation in Sec. IIC. In order to include an often-discussed effect of 
resonant states on an alloy bandstructure Q, one needs to use a more complicate self-energy function, e.g. employing 
the coherent potential approximation, see Ref.9. Last but not least, we note that the potential 5U r in Eq.(2) is 
assumed to be spin-independent, i.e. we have neglected the spin-flip processes induced by the short-range defect. All 
these assumptions restrict the quantitative description of the phenomena resulting from short-range defects therefor 
the results obtained do not related to any concrete case. What this paper does, however, it demonstrates clearly that 
the matrix character of the short-range potential energy in Eq. (3) may change the results essentially. 

Finally, we suppose that the consideration performed will stimulate a reexamination the later obtained results 
both for defects with reduced symmetry in A 3 B 5 semiconductors (see discussion in fL3[ ) and for bandstructure of 
A3_B5-based alloys. Application of such kind of approach for another materials (e.g. A4i?4-based alloys, see [ fl4| ) 
requires a special consideration. In spite of the fact that the obtained results demonstrate some peculiarities of A3B5- 
based alloys, a more detail comparison is restricted due to lack of experimental data for the narrow-gap case under 
consideration (see recent papers |l5| ) . Since the components of the matrix potential can not be fixed in the framework 
of the kp-approach developed, a comparison of such type results with numerical calculations is also of interest. 
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